Looking for a charge asymmetry in cosmic rays 
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Abstract. We combine the data from PAMELA and FERMI-LAT cosmic ray experiments by 
introducing a simple sum rule. This allows to investigate whether the lepton excess observed 
by these experiments is charge symmetric or not. We also show how the data can be used to 
predict the positron fraction at energies yet to be explored by the AMS-02 experiment. 



1. Introduction 

The data collected by PAMELA [lj indicate that there is a positron excess in the cosmic ray 
(CR) energy spectrum above 10 GeV. The rising behavior observed by PAMELA does not fit 
previous estimates of the CR formation and propagation implying the possible existence of a 
direct excess of CR positrons of unknown origins. Interestingly PAMELA'S data show no excess 
in the anti-protons. 

While ATIC [2j and PPB-BETS [3j reported unexpected structure in the all-electron spectrum 
in the range 100 GeV- 1 TeV, the picture has changed with the higher-statistics measurements 
by FERMI-LAT [1] and HESS [5], leading to a possible slight additional unknown component 
in the CR e^- flux over and above the specific Moskalenko and Strong model prediction [61 [7] 
which assumes a single-power-law injection spectrum. 

These interesting features have drawn much attention, and many explanations have been 
proposed: For example, these excesses could be due to an inadequate account of the CR 
background in previous modeling; The presence of new astrophysical sources; They could also 
originate from annihilations and/or decays of dark matter. We refer to [8] for a recent review. 

Whatever the origin of these excesses might be, following Ref. [9], we show how to derive 
constraints about their charge asymmetry. Clearly, this will eventually shed some light on the 
physical nature of their source. 

2. A useful sum rule 

We start by writing the observed flux of electrons and positrons as the sum of two contributions: 
A background component, 4>±, due to known astrophysical sources (which actually are better 
known for the electrons rather than the positrons), and an unknown component, <jr±, in formulae: 

4> ± = 4> U ± + 4>1. (l) 

The component <ft± is the one needed to explain the features in the spectra observed by PAMELA 
and FERMI-LAT. 



These experiments measure respectively the positron fraction and the total electron and 
positron fluxes as a function of the energy E of the detected e^, i.e.: 



p{e) = ,pJ) + Ije) - + #-<*>. W 

The left-hand side of the equations above refer to the experimental measures. The contribution 
from the unknown source is then: 

<j> u + (E) = P{E) F{E) — <j>+{E) , (3) 
4> U (E) = F(E) (1 - P(E)) - <p B (E) . (4) 

In terms of their difference and sum: 

<t>l(E) - <f> u _(E) = F(E) (2P(E) - 1) + (^(E) - <j> B {E)) , 
4>l(E) + <f>V(E) = F(E)-(cp B (E) + ^(E)) . 

(5) 

The latter equation implies F(E) > ^(E) + cj) B (E). We model the background spectrum 
using ^(E) = NbB^(E), where Nb is a normalization coefficient such that F(E)/{B~ (E) + 
B + {E)) > N B and B ± (E) are provided using specific astrophysical models. In this paper we 
adopt the popular Moskalenko and Strong model [7j, for which Nb is less than 0.75 and 
B^(E) are given by: 
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where E is measured in GeV and the Bs in GeV" 1 cm _2 sec _1 sr _1 units. We checked that our 
results remain unchanged when replacing the parameterization above with the one adopted by 
the Fermi Collaboration (model zero) |10| 111]. 

It is convenient to introduce the following parameter: 

/ps _ <fi-(E) F(E) (1 - P(E)) - <£(E) 
ruW <tF + {E) P{E) F(E) - 4>+(E) ' UUj 

which quantifies the level of charge asymmetry in the unknown component of electrons and 
positrons. The equation above can be rewritten as 

Although the sum rule R(E) seems to depend on the energy it should, in fact, be a constant as 
is clear from the right hand side of the previous equation. This leads to a nontrivial constraint 
linking together in an explicit form the experimental results, the model of the backgrounds and 
the dependence on the energy of the unknown components. 



We now turn to the actual data and show in which way the sum rule provides essential 
information on the unknown components of the CRs. 

Since we use simultaneously the results of FERMI-LAT and PAMELA we are obliged to 
consider only the common energy range. Note that the CRs energy range below 10 — 20 GeVs, 
where the spectrum is affected by the Sun, is outside the common range. Within the relevant 
but limited range of energies we will consider here it is therefore sensible to assume rjj to be 
nearly constant. 

We find useful to plot the function R{E) for different values of r\j in order to test the sum 
rule. This would imply that this function is independent of the energy. The associated constant 
value would then be identified with the background CRs normalization factor Nb- We report 
the results in Fig. [TJ The straight (red) line is the Nb = 0.75 value which is the largest one can 
assume for the background not to be larger than the FERMI-LAT results. The shaded region is 
obtained by using the PAMELA and FERMI-LAT data released in 2010, assuming a la error 
band for both. We observe that there is a clear tendency for the combined data to predict a 
lower value of the constant Nb, say about 0.66,0.64,0.62,0.55 for increasing value of the ratio 
rjj = 0, 1, 2,4. This is clear when looking, from top to bottom, at the different panels of Fig. [TJ 
It is interesting to note that we find a plateau, in the relevant energy range, up to rjj near the 
value of 4 when R(E) starts showing some deviation. 




Figure 1. Ratio R{E) as a function of the energy E of electrons and positrons and for values 
of r\j = 0,1,2,4. The shaded region accounts for the la error in PAMELA and FERMI-LAT 
data (2010 release). Secondaries are estimated according to the expressions in [HE]. The value 
of R(E) thus corresponds to the normalization factor for the background, called Nb in the text. 

Given the large uncertainties (for PAMELA, the 2010 la error bands are larger than in 2008) 
we cannot yet provide a more solid conclusion, however we can use the derived normalization Nb 



for each different ratios of the unknown components to predict the positron fraction at energies 
higher than the ones provided so far by PAMELA. 



3. PAMELA and the charge asymmetry of the unknown component 

In order to be able to make such a prediction we first rewrite eq. (jlip as follows: 



P(E) = tt^ y 1 - tw 11 " ru WE) ] ) ' (12) 

where we use for each rjj the estimated associated Nb (referring to Fig. 1, we consider the 
range of Nb values allowed for E ~ 25 GeV). The different predictions for the positron fraction, 
assuming that rjj remains constant over the entire energy range, up to 1000 GeV are shown in 
Fig. [2j The resulting picture disfavours both very small and very large values of t\j which, in 
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Figure 2. Model independent prediction for the positron fraction P(E) as a function of the 
energy E of electrons and positrons and for vales of rjj = 0,1/2,1,2,4, from top to bottom. 
Secondaries are estimated according to the expressions in [61 [7] reported in the main text. The 
shaded regions account for the range of the allowed Nb values (see Fig. 1), obtained in turn by 
considering the 2010 PAMELA and FERMI-LAT data at la. We also display the 2010 (red) 
and 2008 (black) PAMELA data, with la error bars. 

turn, means that one expects the electron fraction to be neither be small nor too large with 
respect to the positron fraction. 

As the experimental data on the positron fraction should hopefully soon become more 
accurate [12], the method outlined above will allow to extract crucial information about the 
charge asymmetry of the unknown component of the electrons and positrons in CRs. 

4. Conclusions 

The model independent analysis of the combined PAMELA and FERMI-LAT data proposed in 
ref. [9] shows that current data still allow for approximately equal contributions of the electrons 



and positrons from the unknown components of the associated CRs, but disfavor electron to 
positron fractions much smaller than 1/2 or larger than 4. Hence, the typical oversimplifying 
model assumption of charge symmetry (rjj = 1) used so far, actually constitutes just a small 
portion of the allowed models still left unconstrained by the present data. 

Interestingly enough, we emphasize that the planned improvements in the experimental 
sensitivity to the positron fraction, as well as further data at energies higher than the one 
explored so far by PAMELA, could reveal whether the unknown component of the CR electrons 
and positrons is charge asymmetric or not. If it turned out to be charge asymmetric, an 
explanation in terms of new physics beyond the Standard Model would clearly become more 
compelling with respect to an astrophysical one. 
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